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By using the method of orthogonal polynomials we ana- 
lyze the statistical properties of complex eigenvalues of ran- 
dom matrices describing a crossover from Hermitian matrices 
characterized by the Wigner- Dyson statistics of real eigenval- 
ues to strongly non-Hermitian ones whose complex eigenval- 
ues were studied by Ginibre. Two-point statistical measures 
(as e.g. spectral form factor, number variance and small dis- 
tance behavior of the nearest neighbor distance distribution 
p{s)) are studied in more detail. In particular, we found that 
the latter function may exhibit unusual behavior p(s) oc s 5 ^ 2 
for some parameter values. 

PACS numbers: 05.45. +b 



Complex eigenvalues of large random matrices have re- 
cently attracted much interest |l|-|lof . Initially, physicists 
used the complex eigenvalues to describe the generic sta- 
tistical properties of resonances in open quantum chaotic 
systems jllJQ and characterize typical features of dissipa- 
tive chaotic quantum maps . Complex and real asym- 
metric matrices were also used to simulate the geomet- 
ric and topological properties of two-dimensional froth 
p3| and to investigate chaotic dynamics of asymmetric 
neural networks I" 



Very recently complex random 
matrices were shown to be relevant for understanding of 
spontaneous chiral symmetry breaking in quantum chro- 
modynamics and also emerged in studies on advection 
of the passive scalar by a quenched random velocity field 
§• 

Traditional mathematical treatment of random matri- 
ces with no symmetry conditions imposed goes back to 
the pioneering work by Ginibre jlTj who determined all 
the correlation functions of the eigenvalues in an ensem- 
ble of complex matrices with independent Gaussian en- 
tries. The progr ess in the field was rather slow but steady, 
see 

Surprisingly, all these studies completely disregarded 
the existence of a nontrivial regime of weak non- 
Hermiticity recognized in our preceding works |l],|5| . This 
regime occurs when the imaginary part of typical eigen- 
values is comparable with the mean separation between 
neighboring eigenvalues along the real axis. In partic- 
ular, the mean density of complex eigenvalues derived 
in H turned out to be described by a formula con- 
taining as two opposite limiting cases both the Wigner 
semicircular density of real eigenvalues of the Hermitian 



matrices and the uniform density of complex eigenval- 
ues predicted for usual non-Hermitian random matrices 
0dT|,0|]. Very recently, Efetov || showed the rel- 
evance of almost-Hcrmitian random matrices to a very 
interesting problem of motion of flux lines in supercon- 
ductors with columnar defects He also managed to 
derive the density of complex eigenvalues for a related, 
but different set of almost-symmetric real random ma- 
trices. This development clearly shows that, apart from 
being a rich and largely unexplored mathematical object, 
almost-Hcrmitian random matrices enjoy direct physical 
applications and deserve to be studied in more detail. 

Although giving an important insight into the problem, 
both papers § and § suffered fr om the same deficiency: 
the present state of art in the application of the super- 
symmetry technique [0] used in gives little hope of 
access to quantities describing correlations between dif- 
ferent eigenvalues in the complex plane due to unsur- 
mountable technical difficulties. At the same time, it is 
just the existence of these universal correlations which 
makes the random matrix theory such a powerful and 
versatile tool of research in different branches of modern 
theoretical physics, see e.g. review p3 |. 

In the present paper we develop the rigorous math- 
ematical theory of almost-Hermitian Gaussian random 
matrices based on the method of orthogonal polynomi- 
als. Such a method is free from the above mentioned 
problem and allows us to study correlation properties of 
the complex spectra to the same degree as is typical for 
earlier studied classes of random matrices. 

To introduce the notion of almost-Hermitian random 
matrices, let us recall that every complex matrix J can 
be decomposed into a sum of its Hermitian and skew- 
Hermitian parts. Following this we consider an ensemble 
of random N x N complex matrices J = H + ivA, where 
H and A are taken independently from a Gaussian Uni- 
tary Ensemble (GUE) of Hermitian matrices with the 

probability density V(X) = Q N X cxp ( - N/2 J§ Tr X 

X = X+. The parameter v, < v 2 < 1, is a natu- 
ral measure of non-Hermiticity in our ensemble. When 
v 2 — 0, J belongs to the GUE and when v 2 = 1, J 
belongs to the ensemble of complex matrices considered 
by Ginibre |T^]. If v 2 is fixed, as N — > oo the eigen- 
values of J are uniformly distributed in an ellipse in the 
complex plane 0,[l4| with the imaginary part being typ- 
ically of the order of Jq. This is just the regime of strong 
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non-Hermiticity. In contrast, the regime of weak non- 
Hermiticity is determined by the condition v 2 = 0(1/ N) 
when N — > oo, see the more detailed discussion in g. In 
other words, we scale the parameter v with the matrix 
dimension N as 2v — a/y/N and consider a to be fixed 
as N — > oo. 

Let us now introduce a new parameter r = (1 — 
v 2 )/(l + v 2 ) and choose the scale constant Jq to be 
equal to (1 + t)/2, for the sake of convenience. The 
parameter r controls the magnitude of correlation be- 
tween Jjk and Jflj-. (JjkJkj) = T /N, hence the degree of 
non-Hermiticity. This is easily seen from the probability 
density function for our ensemble of the random matrices 
J: 



V(J) 



Cjy 1 exp 



N 



Tr( J,P - t Re J 2 )} , (1) 



where C N = [tt 2 (1 - t 2 )/N 2 ] n2 ' 2 . All the J jk have zero 
mean and variance (\Jjk\ 2 ) = 1/N and only Jjk and Jkj 
are pairwise correlated. If r = all the Jjk are mutu- 
ally independent and we have maximum non-Hermiticity. 
When r approaches unity, Jjk and JL are related via 
Jjk = J^j and we are back to an ensemble of Hermitian 
matrices. 

Our first goal is to determine the n-eigenvalue correla- 
tion functions in the ensemble of random matrices spec- 
ified by Eq. (|l|). The density of the joint distribution of 
eigenvalues in the ensemble is given by 



~Pn\Zi, ■ ■ ■ , Zn) 



N N(N+l)/2 



ir N l\---M(l-T 2 ) N / 2 



(2) 



M N 
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One can derive Eq. following more or less standard 
reasoning. By definition, one obtains Vn(Z\, . . . Zn) in- 
tegrating V( J) over the surface of all complex matrices 
whose eigenvalues are Z\ , . . . Zn . We disregard the ma- 
trices with multiple eigenvalues because the set of these 
matrices has zero probability measure of Eq. (^). Ev- 
ery complex matrix with distinct eigenvalues Zi, . . . , Zn 
can be decomposed as J — U(Z + R)U\ where Z = 
diag{Zi, . . . Z^}, U is a unitary matrix, and R is a 
strictly upper-triangular (with zero diagonal) complex 
matrix. In this decomposition Z is determined modulo 
permutation of the eigenvalues and U is determined mod- 
ulo multiplication to the right by a diagonal unitary ma- 
trix. Therefore if we label the eigenvalues and require the 
first non-zero element in each column of U to be positive, 
then the decomposition is unique. We can parameter- 
ize the above mentioned surface by the matrices R and 
U . The Jacobian of the corresponding transformation is 
given by the squared modulus of the Vandcrmonde de- 
terminant [the last factor on the r.h.s. in Eq. @]. The 
integrand V(J) does not depend on U being a Gaussian 
function in R. So, integrating out these two variables is 



straightforward. The resulting expression (with the re- 
quirement of labeling the eigenvalues being removed) is 
Eq. (|). 

The form of the distribution Eq. (EH allows one to cm- 
ploy the powerful method of orthogonal polynomials p8[ . 
Let H n (z) denote the n-th Hermite polynomial, 



rr ( X (±0" ^ 2 



(11 I" exp j — =f izt 



(3) 



The crucial observation borrowed from the paper |23| ] is 
that the polynomials 



Pn(Z) 



r n / 2 VN 



y/^V^.(l - r 2 )V4 IV r 



(4) 



n = 0, 1, 2, . . ., are orthogonal in the complex plane Z 
X + iY with the weight function 



w 2 (Z) = cxpt - 



N 



(1 -T 2 



r {z2 , v .,, 



Z* 



i. e. / d 2 Zp n (Z) Pm (Z*)w 2 (Z) = 5 nm , where d 2 Z = 
dXdY. The standard machinery of the method of or- 
thogonal polynomials |l8| ] yields the n-eigenvalue corre- 
lation functions 



Rn(Z\, Z n 



N\ 



d 2 ! Z, 



n+X- 



d 2 z N r N {z} (5) 



(JV-n)! 
in the form 

R n (Z x , ...,Z n ) = det [K N (Z J ,Zl)] n ] k=1 , 
where the kernel Kjq(Z\ 1 Z%) is given by 



JV-l 



K N (Z 1 ,Z 2 *) = w(Zi)w(Z^) } Pn (Z 1 )p n (Z2) 



(6) 



With Eqs. (|)-@ in hand, let us first examine the 
regime of strong non-Hermiticity, i.e. the case when 
limjv_,.oo(l — r) > 0. In this regime the following Mehler's 
formula wich follows from Eq. (ph is helpful: 



T 

V -rH n (zi)H n (z 2 ) 



: exp 



ziz 2 - - (4 + 



Exploiting this relation one easily shows that in the 
regime of strong non-Hermiticity all the correlation func- 
tions R n (Zi, Z 2 , ■ ■ ■ , Z n ) coincide, after trivial rescaling 
Zj -> (1 - T 2 )- l / 2 Zj, with those found by Ginibre @- 
Now we move on to the regime of weak non- 
Hermiticity: We will demonstrate that in this regime 
new non-trivial correlations occur on the scale: lv&Z\ 2 = 
0(1/ N), ReZi - ReZ 2 = 0(1/ N). Correspondingly, 
we introduce new variables x,yi,y 2 ,ui in such a way 
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that: x = Re(Z 1 + Z 2 )/2, y h2 = Nlm(Z h2 ), uj = 
NRe (Zi — Z 2 ), and consider them finite when perform- 
ing the limit N — > oo. 

Substituting Eq.(^) into Eq.(||) and using the above 
definitions we can explicitly perform the limit N — ► oo, 
taking into account that liniAr^oo N(l — r) = a 2 /2. The 
details of the procedure will be given elsewhere p4|] . As 
a result, the kernel Kn(Z\, Zj*) is given by the following 
expression: 



K N {Z x ,Z*)=g a {y-iuj/2) x 
ix(yi - y 2 



N / y 2 + 

— exp <^ j 

ira or 



+ 



where we have intoduced introduced the notation y 
(yi +2/2)/2 and 



with v sc (x) 



27T 



exp 



2uy 



(8) 



standing for the Wigner semi- 



circular density of real eigenvalues of the Hermitian part 
H of the matrices J. 

Equation (Q) constitutes the most important result of 
the present publication. The kernel Kn given by Eq. (|t]) 
determines all the properties of complex eigenvalues in 
the regime of weak non-Hermiticity. For instance, the 
mean value of the density p{Z) = J2i=i S {2) (Z - Zi) of 
complex eigenvalues Z = X + iY is given by (p{Z)) = 
K N {Z, Z*) Putting yi = y 2 and u = in Eqs. (@)-(f|) we 
immediately recover the density of complex eigenvalues 
of almost-Hermitian random matrices obtained for the 
first time in our preceding publication Q and reproduced 
recently in |J. 

One of the most informative statistical measures of the 
spectral correlations is the 'connected' part of the two- 
point correlation function of eigenvalue densities: 

(p(Z 1 )p(Z 2 )) c = {p{Z x )) S^(Z 1 - Z 2 ) - Y 2 (Z 1 ,Z 2 ), (9) 

with the so-called cluster function Y 2 (Z\, Z 2 ) being equal 
to Y 2 (Z U Z 2 ) = \K N (Z u Zl)\ 2 . In particular, it deter- 
mines the variance £ 2 (£>) = (n(D) 2 ) — (n(D)) 2 of the 
number n = J D d 2 Z p(Z) of complex eigenvalues in any 
domain D in the complex plane. It also allows one to 
find the small-distance behaviour of the so-called nearest- 
neighbour distance distribution : p(zq,s) , i.e. the prob- 
ability density to have an eigenvalue at a point z of the 
complex plane such that its closest neighbour is at the 
poi nt z\ separated by the distance \z\ — z§\ = s/N, see 

Ml§- 

We see, that in the regime of weak non-Hermiticity the 
cluster function depends on the variable x = Re(Zi + 
Z 2 )/2 only via the semicircular density v sc {x). The pa- 
rameter a = 2nv sc a controls the deviation from Her- 
miticity. When a — > the cluster function tends to GUE 



form Y 2 (yi, 2/2, uj) = N 2 8(y 1 )S(y 2 ) 



z (x)u>) / {-Ku>y 



whereas in the opposite case a» 1 the limits of integra- 
tion in Eq.(||) can be effectively put to ±00, and the cor- 
responding Gaussian integration yields the Ginibre |T^ ] 
result: Y 2 (Z 1 ,Z 2 ) = {N/ira) 2 exp{-iV 2 |Zi - Z 2 \ 2 /a 2 }. 

Fourier-transforming the cluster function over all ar- 
guments yi , y 2 and u> we find the following expression for 
the spectral form-factor. 

2 

b( gi ,q 2 ,k) = N 2 exp{-^- + q2 + 2fc2 ) } X (10) 



TTa 2 (q 1 +q 2 ) 



1M 

2tt 



6{2kv sc (x) - |fc|) 



na 2 (qi + q 2 )/2 

where 9{u) — 1 for u > and zero otherwise. 

As we already mentioned, the knowledge of the cluster 
function allows one to determine the variance of the num- 
ber of eigenvalues in any domain D of the complex plane. 
In the general case this expression is not very transpar- 
ent, however. For this reason we present the explicit form 
for the simplest case, choosing the domain D to be the 
infinite strip of width L (in units of mean spacing along 
the real axis A = {vsc^N) -1 ) oriented perpendicular 
to the real axis: < ReZ < LA; —00 < ImZ < 00. 
One finds: 



E 2 (L) = L 



1 dk 
k 2 



(1 - k) sin 2 (irkL)e 



_ (afi)2 



(11) 



As long as a <C L the function £2 (L) slightly deviates 
from the corresponding expression for the Hermitian ma- 
Let us briefly discuss the gross features of 
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trices 

the number variance at a > 1 when we expect essential 
differences from the Hermitian case. Normally one is in- 
terested in large L asymptotics: L> 1. In that case the 
upper limit of the integral in Eq.([ll|) can be set to infin- 
ity. This implies that in a large domain 1 <C L ~ a the 
number variance grows linearly like = Lf(2L/a) 

with 



dte-*' 



For u — 2L/a <C 1 we have simply / w 1. For u ^> 
1 we have / ~ (y/iru)^ 1 . Thus, £ 2 (L) slows down: 
Tj 2 (L) ss a/2y/ir. Only for exponentially large L such 
that ln(L/a) > a does the term proportional to k^ 1 in 



the integrand of Eq. (jllj) contribute significantly. Then it 
leads to the slower logarithmic increase: 



S 2 (£) 



2y / 7T 7T 



no 



corrections, 



where 7 is Euler's constant. This behaviour is reminis- 
cent of that typical for real eigenvalues of the Hermitian 
matrices. 

One can also derive the following expression for small- 
distance behaviour of the probability density to have one 
eigenvalue at the point zq = X + iyo/N and its closest 
neighbor at the distance \z\ — Zq\ = s/N: 



3 



2s 3 / dBexp 



■9a(2/o)^lnffa(yo) _ 2 v| 

' ! = (Wp 6 " 



(12) 



: (y - s cos ( 



It is easy to see that performing the limit a — ► 
we are back to familiar GUE quadratic level repulsion: 
p Q ^o(^o, s < 1) oc S(yo)s 2 1 whereas in the opposite limit 
a — > oo one has: p a ^oo(zo, s <C 1) = ^(s/a 2 ) 3 in agree- 
ment with the cubic repulsion generic for strongly non- 
Hermitian random matrices |p7y^2| j7[|. In general, the ex- 
pression in Eq.(12) describes a smooth crossover between 
the two regimes. 

To this end, an interesting regime occurs for the 'ob- 
servation points' zq situated sufficiently far from the real 
axis: \y a \ 3> 2~ 3 / 2 a. It turns out that for such points the 
distribution p(zo,s) displays an unusual s 5 / 2 power-law 
behaviour: 



P a (z ,s < 1) 



gl (vo)-§^ in g a {yo ) 

2(2^)3/2a| 2/0 | 1 /2 



- 4 3 s 5 / 2 



(13) 



in the parametrically large domain of distances: 
a 2 /4 1 y 1 <C s -C 2\y \. Unfortunately, the low density 
of complex eigenvalues in such points (reflected by the 
presence of the Gaussian factor in the expression above) 
would make any direct verification of such an unconven- 
tional behaviour to be a difficult numerical problem. 

Finally, we would like to make the following remark. 
So far all the expressions were derived rigorously for the 
Gaussian almost-Hcrmitian random matrices by using 
the probability measure of Eq.(^) as a starting point. 
However, employing the supersymmctry technique one 
is able to give strong, albeit not mathematically rigor- 
ous arguments, [^4| that all the results obtained in this 
way are fairly universal and equally applicable for a much 
broader class of almost-Hermitian matrices J = H + ivA, 
provided the semicircular density v sc (x) is replaced by 
the actual eigenvalue density v(x) of the Hermitian ma- 
trix H on the real line at point x. 

In conclusion, we have developed a rigorous mathemat- 
ical theory of almost-Hermitian Gaussian random ma- 
trices based on the method of orthogonal polynomials. 
This has enabled us to study in great detail the two- 
point correlations between complex eigenvalues in the 
novel crossover regime between Hermitian and strongly 
non-Hcrmitian random matrices. 
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